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Strong spatial dispersion in wire media in the very large wavelength limit 
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It is found that there exist composite media that exhibit strong spatial dispersion even in the 
very large wavelength limit. This follows from the study of lattices of ideally conducting parallel 
thin wires (wire media). In fact, our analysis reveals that the description of this medium by means 
of a local dispersive uniaxial dielectric tensor is not complete, leading to unphysical results for 
the propagation of electromagnetic waves at any frequencies. Since non-local constitutive relations 
have been usually considered in the past as a second order approximation, meaningful in the short 
wavelength limit, the aforementioned result presents a relevant theoretical interest. In addition, 
since such wire media have been recently used as a constituent of some discrete artificial media (or 
metamaterials) , the reported results open the question of the relevance of the spatial dispersion in 
the characterization of these artificial media. 
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Causality imposes that all material media must be dis- 
persive. In most cases this behavior results in local dis- 
persive constitutive relations - i.e. in frequency depen- 
dent constitutive permittivity and permeability tensors. 
Non-local dispersive behavior (i.e. spatial dispersion), 
which results in constitutive operators depending also on 
the spatial derivatives of the mean fields (or, for plane 
electromagnetic waves, on the wavevector components) 
is usually considered as a small effect, meaningful in the 
short wavelength limit. Specifically, spatial dispersion 
will always appear when the higher order terms in the 
series expansion of the constitutive parameters in power 
series of the dimensionless parameter a/X (a is the lattice 
constant of the crystal and A the wavelength inside the 
medium) are not neglected jJ]. Thus, it is usually as- 
sumed that non-local dispersive reations are only mean- 
ingful when A approaches to a. The usually weak natural 
optical activity of some materials is a well known exam- 
ple of the application of this principle [pj. When such 
principle is translated to the analysis of discrete artificial 
media, also called metamaterials, it would imply that 
non-local dispersive constitutive relations are only ex- 
pected to be a small refinement of the local constitutive 
relations usually considered. However, there is at least a 
counter example for this assumption. 

The parallel wire medium is a medium formed by a reg- 
ular lattice of ideally conducting wires with small radii 
compared to the lattice periods and the wavelength, see 
Fig. |} It has been known in microwave applications for 
a long time ^ [| |j as an artificial dielectric, also called 
rodded medium and quasistatic models of the effective 
permittivity are available [^). Recently, some attention 
to wire media has been paid also in optics (e.g. [|[ [?]]) 
and in the realization of left-handed media) [pi, M as com- 



posite media made from lattices of long conducting wires 
and split ring resonators j^oj [Til , (see discussions 

corresponding to that in |Q\3|| ) . The electromagnetic re- 
ponse of the specific wire medium shown in Fig. ^ is 
analyzed in || fj]] following different approaches. Both 
analysis are carried out for wave propagation perpendic- 
ular to the wires and show that, for electric field polar- 
ization along the wires, the medium is characterized (if 
a/A , b/X <C 1) by a frequency dependent effective dielec- 
tric constant given by: 

The constant too (the corresponding wavenumber ko — 
LUoy/^ofJ-o) i n epl plays the role of an equivalent "plasma 
frequency." Thus, this medium is often called "artificial 
plasma" since the ideal (collisionless) electron plasma is 
described by the same equivalent parameter. If the wires 
are assumed to be very thin, so that their polarization 
in the direction orthogonal to the wires can be neglected, 
the effective permittivity for electric field polarization or- 
thogonal to the wires is £o- 

The aforementioned analysis suggest that this wire 
medium could be modeled as an uniaxial dielectric with 
the following local permittivity dyadic: 

I = ez z +£o(x x + y yo), (2) 

whose permittivity in the axial direction, e, would be 
given by epl. However, it will be shown in the following 
that this naive hypothesis lead to unphysical results and 
must be substituted by a non-local dispersive relation. 
In fact, assuming that the medium can be described by 
the uniaxial dyadic uni, the dispersion equation for ex- 
traordinary plane waves (E z ^ 0) with the wavevector 
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FIG. 1: The geometry of wire media: a lattice of parallel 
ideally conducting thin wires. 



(q x , q y ,q z ) T in this uniaxial dielectric reads jT^, |l5[ ]: 



eo(9x + 9«) = £ ( fc2 -^), 



where k — aj^/eojlo is the phase constant of the host ma- 
trix. On the other hand, these extraordinary waves corre- 
spond to the well known TAf-(to z) set of modes, allowed 
by the invariance of the boundary conditions along z. 
Thus, for any extraordinary wave travelling with a phase 
constant q z along the z-axis, the E z field must satisfy the 
Helmholtz equation 



d | d | 2 
dx 2 dy 2 



qt)}E z =0, 



rationale suggests that the considered wire media still 
can be described by the permittivity dyadic uni, but the 
axial permittivity e must be a non-local parameter of the 
form eps. The conventional expression epl would be only 
a particular case of eps, valid for wave propagation in the 
x — y plane. 

The main difference between the local uniaxial model, 
epl, and the proposed non-local model, eps, for the par- 
allel wire medium is that the non-local model predicts 
a stop band (at frequencies below luq = fco/-i/EoMo) f° r 
extraordinary waves propagating along any direction in 
the media. On the contrary, epl-disuni predict propa- 
gation of extraordinary waves at any frequency provided 
q z > k = uj^/HqEq, Thus, both models predict qualitative 
very different behaviors, even near the cutoff "plasma" 
frequency, cj , where q 2 — > (i.e. a/A — ► 0). That is, the 
non-locality of the proposed constitutive relations affects 
the electromagnetic response of the medium even in the 
very large wavelength limit, thus being important for any 
values of the a/ A ratio inside the medium. Other rele- 
vant differences between the predictions of both models 
will be developed along this letter. 

The rigorous proof of eps is based on the local field 
approach which is described in detail in jl6|. In Jll| the 
low frequency stop band of the wire medium has been 
analyzed, as well as its high frequency band gap struc- 
ture. This analysis reveals that, in the thin wire medium 
(Fig. |l|) and for q z ^ k, two sets of modes can propa- 
gate: ordinary (with E z = 0) and extraordinary (with 
E z ^ 0) waves. The ordinary waves do not interact with 
the wires and propagate in the host media. For extraor- 
dinary waves, an explicit dispersion equation connecting 
the wave vector q = (q x , q y , q z ) T with the wave number 
of the host isotropic matrix k = aj^/£oMo bas been derived 
( 4 ) in (H). It can be written as: 



(3) 



with the boundary condition E z = on the wires. It is 
clear from this equation that any "plane" extraordinary 
wave must satisfy 



Klx,q y ,q z ) = \Jk 2 {q Xl q Vl 0) + q 2 



(5) 



This last result is incompatible with epl-disuni, as can be 
easily seen by substitution of epl into disuni. However, if 
we choose 



e(k,q z 
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k 2 



k 2 - q\ 



(6) 



instead of epl, then uni and disuni become compatible 
with total, giving the following dispersion equation for 
the plane wave: 

q 2 = q 2 x + q 2 y + q 2 z = k 2 -k 2 , (7) 

where we have assumed that q z ^ k (the case with q z — k 
will be analyzed at the end of this letter). The above 
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Here k x denotes the x-component of n-th Floquet mode 
wave vector: 



fe(") = -jy^ + ^j +q 2 -k 2 , Re{v^)}>0. 

(9) 

The other notations are clear from Fig. |l|. Numerical so- 
lution of this dispersion equation shows that there exists 
a low frequency stop band for all propagation directions 
(except for the particular case of q z = k that will be 
analyzed latter). Let us simplify the dispersion equa- 
tion disp for the quasi-static case a, b <C w/k. Using the 
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Taylor expansion of sin(a;) and cos(cc) functions for small 
arguments we obtain qq with 
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where s = Vab, r = a/b and 
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Therefore, we have shown the suitability of the suggested 
approach for the description of the wire medium, with k 
given by kO. Parameter kg here plays the role of the wave 
number corresponding to the plasma frequency. More 
exactly, it indicates the upper edge of the low frequency 
stop band. Naturally, fco as a function of two lattice pe- 
riods a and b is a symmetric function: ko(a, b) = ko(b, a). 
It means that function F(r) has the following property: 
F(l/r) = F(r). For the commonly used case of the 
square grid (a = b), F(l) — 0.5275. Expression kO looks 
similar to the approximate expressions for the plasma fre- 
quency developed earlier in |2| ||, ||, ^| , but for thin wires 
it is more accurate and takes into account the geometry 
of the lattice. Notice that the dispersion equation qq for 
extraordinary waves is indifferent to the direction of the 
wires axis: the wave vector components q x -,qy^qz enter 
into this equation completely symmetrically. It means 
that, within the low frequency stop band, the extraor- 
dinary wave decays with the same decay factor along all 
directions in space. The same can be said for propagation 
in the first frequency pass band. This isotropy of the dis- 
persion equation is rather surprising since the medium is 
strongly anisotropic. However, it can be shown from the 
very fundamental facts summarized in total and disuni 
by assuming, as usual, that epl is valid for extraordinary 
waves propagating in the x — y plane. 

It is possible to transit eps from the spectral domain 
(q,w) to the physical domain (r,t). The following non- 
local material equation can be derived from eps using the 
double Fourier transform: 
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where c = 1/ ^/£oMo is the speed of light in the host ma- 
trix. Here, the area of integration in the z — t plane is 
the light cone \z — z'\ < c(t — t'). In other words, the 
kernel in the Fourier convolution is u[c(t — t') — \z — z'\], 
where u(x) is the Heaviside step function. It means that 
the point (a;, y, z) inside the wire medium (described as 
a dispersive continuum) at moment t is affected by the 



^-components of electric fields coming from the domain 
(x, y,z± c(t — i')) surrounding (along the wire axis) this 
point during all the past time (f < t). This result il- 
lustrates the consistency of the reported model from the 
relativistic standpoint. 

In the following we will describe some relevant effects 
in the analyzed parallel wire medium, associated with 
the spatial dispersion. Refraction and reflection of plane 
waves at a plane interface show strong differences be- 
tween the local and non-local models. Let us consider 
an interface between an isotropic dielectric with the per- 
mittivity £j and a uniaxial dielectric with e described 
by epl. The interface is in the (y — z)-plane and it is 
illuminated by a plane wave coming from the isotropic 
dielectric. The wave vector and electric field vector lie in 
(x — z) plane (q y = 0, E y = 0). The incidence angle of the 
plane wave is 9. If e\ > Eq, e < 0, and sin 2 (#) < eq/ex the 
wave will be completely reflected, but for sin (6) > eq/ex 
some part of the wave will be transmitted through the 
interface. This transmitted wave will be an extraordi- 
nary wave, as it follows from its polarization state, and 
can be excited at any frequency. This amazing behavior 
disappears when the non-local model summarized in eps 
is used. Indeed, if the non-local axial permittivity eps 
is used for the wire medium, we observe that no trans- 
mission inside the wire medium is possible for k < feo- 
At k — ko transmission is possible only in the case of 
the normal incidence. Only if k > fco, a refracted wave 
appears. 

Let us next consider the guidance of electromagnetic 
waves in a parallel-plate waveguide infinite in the x and 
y directions and bounded by parallel perfectly conduct- 
ing planes orthogonal to the z-axis. Separation between 
the conducting walls is d. We assume that this waveg- 
uide is filled with a wire medium with the wires along 
the z-direction. We will consider eigenwave propagation 
along the x axis of the TMoi mode {H yi E Xl E z ^ 0). 
For waveguides filled by a local uniaxial dielectric with 
anisotropy axis along the z direction we have from disuni 



£ ql = e(k 2 - q 
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(13) 



If e > Equation propuniy gives a cut-off for k < ir/d 
and propagation for k > ir/d. In contrast, if e < 0, 
propagation is allowed when k < ir/d (and forbidden 
for k > n/d). Within this pass band a backward wave 
(dq/duj < 0) propagates, as one can see in Fig. [2] (thin 
solid lines). 

This amazing behavior disappears if one fills the 
waveguide with the analyzed non-local wire medium. Us- 
ing qq, we have in this case 



t + ll = k2 - k l q* = \lk 2 -k 2 



(14) 
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and we obtain the usual frequency behavior: cut- 
off for k < yj (ir/d) 2 + k 2 and propagation for k > 
\J (ir/d) 2 + fcg. An increase of the cut-off frequency is 
observed compared to the case when there is no filling 
medium, as one can see in Fig. || (thick solid line and 
dashed line). 




Normalized frequency, klk^ 

FIG. 2: Normalized propagation factors in a parallel-plate 
waveguide vs. normalized frequency for different types of 
waveguide filling: thin solid lines - uniaxial dielectric with a 
negative permittivity; thick solid line - wire medium; dashed 
line - empty waveguide. Fhe wire medium and the uniaxial 
dielectric have the same ko = n/(2d). 

Let us finally analyze the propagation of plane waves 
along the wire medium for the particular case of q z = k. 
In this case, the non-local permittivity along the 2-axis 
eps becomes infinite. To avoid the singularity problem 
we use material equation of the form E = e 1 D. In 
this case the Maxwell equations have plane wave solu- 
tions for all frequencies. For these waves the transverse 
wave-vector = (q x ,Qy) T is arbitrary. The waves are 
transverse with respect to the wire axis: H z = and 
E z = 0. The electric field is parallel to the transverse 
wave vector, x E = 0. 

Such waves can be interpreted as transmission-line 
modes propagating along the parallel wires. In fact, a 
set of N infinite parallel wires can be viewed as a system 
of coupled transmission lines. This system can support N 
degenerate transmission line waves with H z = 0, E z = 
and phase constant q z = k. The electric field of these 
waves can be obtained from 

E= -(u x d x + u y dy)<l>(x,y)exp(-jkz), (15) 

where 4>(x, y) is a quasi-electrostatic potential taking con- 
stant but arbitrary values <j) n (n = 1,2, N) at each 
wire. In fact, the plane wave with transverse wave- 



number q^ and q z = k corresponds to the transmission- 
line wave with <p n = <po exp (— jqj_ • r n ), where r n = 
(x n ,yn) T is the location of the n-th wire in the trans- 
verse (x — y) plane. 

In summary, it has been shown that parallel wire media 
possess very strong spatial dispersion effects at any fre- 
quencies, including the very large wavelength limit. An 
analytical model for the non-local permittivity dyadic of 
these media has been presented and discussed. Inconsis- 
tency of the local model for parallel wire media with non- 
vanishing wave-vector component along the wires has 
been shown. Dramatic differences in the predicted behav- 
ior of that media, arising from the use of the conventional 
local and/or the non-local model for the permittivity are 
shown. Finally, the proposed non-local model for the per- 
mittivity has been found to be also suitable for the de- 
scription of the transmission-line modes of the structure. 
We feel that the reported results opens the question of 
the role of spatial dispersion in the adequate characteriza- 
tion of discrete metamaterials as effective media, at least 
if arbitrary directions of propagation and/or polarization 
of the electromagnetic field should be considered in the 
analysis. In addition, an example has been presented of 
an effective medium in which spatial dispersion is impor- 
tant at any frequency, in contrast with some commonly 
assumed ideas about the physical relevance of this effect. 
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